Transition to amplitude death in scale-free networks 
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Transition to amplitude death in scale-free networks of nonlinear oscillators is investigated. As the coupling 
strength increases, the network will undergo three stages in approaching to the state of complete amplitude 
death. The first stage is featured by a "stair-like" distribution of the node amplitude, and the transition is 
accomplished by a hierarchical death of the amplitude stairs. The second and third stages are characterized 
by, respectively, a continuing elimination of the synchronous clusters and a fast death of the non- synchronized 
nodes. 
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The discoveries of the small-world and scale-free properties 
in natural and man-made systems have brought a new surge to 
the studyof collective behaviors in coupled dynamical sys- 
tems mM [U, wnere th e scope of the traditional studies had 
been significantly extended and a number of new phenomena 
had been identified ||4j] . A typical example could be the syn- 
chronization of populations of coupled nonlinear oscillators. 
Recent studies have shown that, due to the shorted network di- 
ameter, small-world networks in general have the higher syn- 
chronizability than regular networks of the same parameters 
lUt]. In studying scale-free networks, it has been found that 
synchronization is much influenced by the few large-degree 
nodes than the majority small-degree ones [6]. Nowadays, 
the interplay between the network topology and dynamics had 
been one of the focusing issues in nonlinear studies, and there 
are many questions waiting for explorations ]4|]. 

Amplitude death (AD) refers to the cessation of oscilla- 
tion of coupled oscillators when their parameters are consider- 
ablymismatched or there exists time delay in their couplings 
I0M Htl- Typical examples include the arrhythmia of the 
cardiac pacemakers fllOll. the quenched oscillation of coupled 
electronic circuits fllll . optical oscillators IU2I1 . and chemical 
reactors 111 311 . Recently, this study has been extended to spatio- 
temporal systems, where the influences of the frequency dis- 
tribution and network topology have been addressed. A chain 
of nonlinear oscillators of monotonic frequency distribution 
has been studied in Ref. ltl4ll . where partial amplitude death 
(PAD) (a state where only partial of the network nodes are 
dead) has been observed before the complete amplitude death 
(CAD) of the system. The authors have also found that, by 
disordering the frequency distribution, the generation of CAD 
can be significantly postponed. Eliminating CAD by disor- 
dered network structure has been studied in Ref. [16], where 
the small-world networks have been found to be more resis- 
tive to AD than the regular and random networks. The route 



to CAD in an array of oscillators has been investigated in 
JT5I1 . where the transition is found to be divided into differ- 
ent stages. 

The papers cited above, however, deal with only the case 
of homogeneous networks, i.e., nodes in a network have the 
similar degree. As practical systems usually take the form 
of scale-free networks characterized by the existence of few 
very large-degree nodes, it is intriguing to see how the hetero- 
geneous distribution of the node degree will affect the genera- 
tion of AD. In the present work, by investigating the transition 
of scale-free networks to CAD, we shall address the important 
role of node degree played in AD generation. 

We consider N coupled Landau-Stuart oscillators of the 
following form: 
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Zi(t) = [7 
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where i — 1, 2, . . . , TV is the node index. Zi(t) is a complex 
number denoting the state of the ith oscillator at time t, e is 
the uniform coupling strength, u>i is the natural frequency of 
node i, and r is the oscillation growth rate. The connections 
of the oscillators are defined by matrix A = {fly}, where 
an = 1 if i and j are connected, and ay- 



degree of the ith node is fc, = 



N 
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= otherwise. The 
! <iij. To differentiate the 
node dynamics, uji is randomly chosen from range [wi,^]. 
Without coupling, the trajectory of each oscillator will settle 
to a limit circle of radium | Zi (i) | = r. In our study, we will fix 
all other parameters in Eq. ([TJi, while increasing e to realize 
the transition. 

To measure the degree of the network death, the following 
two macroscopic quantities will be employed: The normal- 
ized network "incoherent" energy, E = (X)j=i \Zi{t)\ 2 ) /N , 
and the number of the "dead" nodes Nd- Here (•) denotes the 
time average over time period T. It is important to note that, 
due to the noise perturbations, in practice the averaged ampli- 
tude of an oscillator, Zi = (\Zi(t)\), can not be exactly zero, 
even through the condition of AD is fulfilled. So it is neces- 
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FIG. 1 : (Color online) The transition to CAD for a scale-free network 
of N — 1000 Landau-Stuart oscillators, (a) The variation of E as 
a function of e, where the transition is roughly divided into three 
stages by the critical couplings e' ~ 0.16 and s" ~ 0.2. The fitted 
exponents of the first and third stages are —9 ± 0.1 and —335 ± 1, 
respectively, (b) The variation of Nd as a function of e. Two different 
amplitude thresholds, Zo = 1 x 10~ 2 (the red curve) and Zo = 
1 x 10~ 10 (the black curve), are used separately to characterize the 
transition. 



FIG. 2: (Color online) For the same network used in Fig. 1, typical 
amplitude distributions observed in the first stage of the transition, 
(a) The distributions are generated by coupling strengths, from top 
to bottom, e = 0.05,0.09,0.1,0.15,0.16. Inset: The inverse re- 
lationship between the critical coupling Ek and the node degree k. 
The symbols are the numerical results, and the solid line is drawn by 
the theory, (b) For the same set of couplings as in (a), the averaged 
stair amplitude (Zj,) versus the node degree k. Each data of (b) is 
averaged over 50 network realizations. 



sary to define some threshold, Zo = (\Z\) <C 1, and regarding 
nodes of Zi < Zq as dead. Apparently, a state of smaller E 
and larger Nd corresponding to a higher degree of amplitude 
death. 

We first study the transition by means of the macroscopic 
quantities E and Nd- The network is generated by the stan- 
dard BA model [2], which consists of N — 1000 nodes and 
has average degree (k) = 6. The degree distribution follows 
roughly the power-law scaling P(k) ~ k~ 3 . To facilitate the 
analysis, we reorder the network nodes by an ascending order 
of their degrees. Thus k\ = k m i n and k^ = k max . For the 
local dynamics, we set r = 0.5 and randomly choose w,; from 
the range [1, 50]. To simulate, we initialize the network with 
random conditions and evolve it according to Eq. (|T). After a 
transient period T' = 1 x 10 3 , we start to calculate E and Nd, 
which are averaged over another period T = 500. The varia- 
tion of E as a function of e is plotted in Fig. 1(a). It is seen that 
the evolution of E can be roughly divided into three stages. 
In the range e £ [0, 0.16], E is decreased exponentially, but 
with a relatively slow speed; in the range e € (0.16,0.2), 
the decrease of E is fasted, but the dependence of E on e is 
complicated (as a matter of fact, the value of E is sensitively 
dependent on the frequency distribution of the oscillators); in 
the range e £ [0.2, 0.25], E is decreased exponentially again, 



but with a much faster speed in comparison with the previous 
stages. The variation of Nd as a function of e is plotted in Fig. 
1(b). For a small amplitude threshold, e.g. Zo = 1 x 10 -10 , 
Nd is found to be smoothly increased from to N within a 
narrow range of e. However, if the amplitude threshold is not 
too small, e.g. Zq > 1 x 10~ 2 , a non-smooth stair-like evolu- 
tion of Nd will be presented. 

To explore the staged transition and the stair structures ap- 
peared in Fig. 1, we proceed to investigate the system dy- 
namics at the microscopic level. Specifically, we shall calcu- 
late the distribution of the node amplitude, {Zi}, and see how 
this distribution is evolved during the transition. In the first 
stage, the node amplitude is evolved as follows. At very small 
coupling e ks 0, all nodes own the same amplitude Zj = 7-, 
regardless of the difference of the node degree. Then, as s 
increases, all the amplitude will be decreased, but with very 
different speed. Specifically, the decrease of the amplitude of 
a node is proportional to the node degree. So the amplitude 
distribution is gradually curved and replaced by a non-smooth 
"stair-like" distribution, say, for example, the one generated 
by e = 0.05 in Fig. 2(a). In this stair-like distribution, nodes 
of the same degree will present the same amplitude, except 
some rare "bursts" which are caused by synchronous clusters 
(to be explained later). Because of the heterogenous degree 
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distribution, the stairs are of different length. More interest- 
ingly, the height of the stairs is gradually stepping down with 
the node degree. That is, the highest stair consists of only the 
smallest-degree nodes and the lowest stair consists of only the 
largest-degree nodes. Increasing e further, the amplitude of 
the largest-degree nodes is hardly decreased, while the ampli- 
tude of the other nodes will be decreased continuously. This 
results a hierarchical death of the amplitude stairs. The pro- 
cess of stair elimination is shown in Fig. 2(a), where some 
typical amplitude distributions observed in the first stage are 
presented. Finally, at about the coupling strength e', the last 
stair, i.e. the stair consisting of the smallest-degree nodes, will 
be eliminated, and the first-stage transition is completed. For 
instance, the distribution generated by s = 0.16 in Fig. 2(a). 
So the stair structures appeared in Fig. 1(b) is just a reflection 
of the hierarchical death of the amplitude stairs. 

The formation of amplitude stairs could be analyzed by 
the mean-field approximation. Noticing that Eq. (Q]l can be 
rewritten as Z, t (t) = [(r 2 - ski) + icjj - \Zi\ 2 ]Zi(t) + G t , 
where G t = e J2i= 

the collective coupling received 
by node i. When e is small, most oscillators of the net- 
work will behave incoherently, therefore G L is small and neg- 
ligible to the node dynamics, especially for the large-degree 
nodes. With this concern, the node dynamics can be simpli- 
fied to Zi(t) = [{r 2 - ski) + i^i - \Z t \ 2 ]Z t (t). By requiring 
\Zi(t)\ = 0, we obtain 

\Z t \ 2 = r 2 -sk l . (2) 

Eq. (0 tells that, for a given coupling strength, the amplitude 
of a node is only dependent on the node degree. This explains 
why the height of the stairs in Fig. 2(a) is decreased with the 
node degree. From Eq. (fJI we can also estimate the averaged 
amplitude of each stair Z}., which is verified by the numerical 
data of Fig. 2(a). Moreover, by setting \Zj\ = OinEq. (0, we 
can obtain the critical coupling Ek where the /c-degree stair is 
eliminated: Ek ~ r 2 /k. The inverse relationship between Zf~ 
and k is verified in the inset plot of Fig. 2(a). 

It should be mentioned that the above analysis works for 
only the case of weakly coupled networks. If the coupling 
strength is not too weak, some nodes in the network could 
behave coherently. In such a case, the value of Gi will be de- 
viated from 0, and the flat stairs predicted by Eq. (fJJ will be 
disturbed. This is just we have observed in simulations. In 
Fig. 2(a), despite the changes of the coupling strength, there 
always exist some bursts in the distributions. These bursts, 
which have random locations and various amplitude, are di- 
rectly resulted from the synchronized nodes. The picture is 
the following. Since the node frequency is randomly chosen, 
there could be the situation that some connected nodes in the 
network have very small frequency mismatch. As the cou- 
pling strength increases, these nodes will be easily synchro- 
nized and form some synchronous clusters (phase synchro- 
nization ifnll ). Once synchronized, nodes will be efficiently 
protected from AD Il4,[l5ll . showing as the amplitude bursts. 
The synchronized nodes, however, may have different ampli- 
tude, as they could be embraced by different set of neighbors. 
To smooth the amplitude bursts caused by synchronization, 
we have calculated the averaged stair amplitude, (Z).), as a 
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FIG. 3: (Color online) For a slightly modified network model, the 
amplitude distributions observed in the second and third stages of 
the transition, (a) e = 0.18 and (b) e = 0.19 are within the second 
stage, (c) e = 0.21 and (d) e = 0.22 are within the third stage, 
(e) For the above amplitude distributions, the variations of (Z^) as a 
function of k. The missed data in the last distribution is caused by 
the limited computing precision. 



function of k. Now a smooth, hierarchical transition of AD is 
presented [Fig. 2(b)]. 

We go on to investigate the transition in the second stage. 
At the end of the first stage, a number of synchronous clus- 
ters are formed in the network. Because of synchronization, 
the amplitude of the synchronized nodes are apparently larger 
than the non-synchronized ones. In the second stage, as the 
coupling strength increases, both the number and the size of 
the synchronous clusters are decreased. Interestingly, it is 
found that the robustness of a cluster is mainly determined 
by the frequency mismatch between the synchronized nodes, 
while is less affected by the node degree or cluster size. To 
show this point more clearly, we have slightly modified the 
network by artificially adjusting the frequency mismatch be- 
tween two connected nodes, (147, 232), to be Suj = 1 x 1CP 3 . 
This pair of nodes, both having the smallest degree k m i n = 3, 
are synchronized at a very small coupling and are kept at large 
amplitude till the very end of the transition. As shown in Fig. 
3, despite of the increase of the coupling strength, the ampli- 
tude of these two nodes is always apparently larger than the 
resting nodes. Different to the situation of the first stage, in 
the second stage the node amplitude is less dependent of the 
node degree, as indicated by Fig. 3(e). 

The second-stage transition is ended up with a state of very 
few synchronous clusters. Typically, each cluster consists of 
only several nodes which have very close natural frequencies, 
like the pair of nodes discussed above. The few survival clus- 
ters, however, are extremely robust and could stand for a very 
large coupling strength. The robust clusters result the new 
form of transition in the third stage: A fast death of the non- 
synchronized nodes accompanied by a slow death of the syn- 
chronized nodes. This property of network transition can be 
read from Fig. 3(b) and (c), where the amplitude of nodes 147 
and 232 are kept at large values, while the amplitude of the 
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FIG. 4: Typical amplitude distributions (R4 = (y/xffij + yf(t))) 
observed in the transition of the network of chaotic Rossler oscilla- 
tors, (a) e = 0.039 and (b) e = 0.05 are within the first stage; (c) 
e = 0.055 and (d) e = 0.06 are within the second stage; and (e) 
e = 0.065 and (f) e = 0.07 are within the third stage. 



other nodes is extremely small (see also Fig. 3(e)). Like the 
second stage, in the third stage the node amplitude is also less 
dependent on the node degree, as can be found from Fig. 3(e). 
With the elimination of the most robust cluster, the whole tran- 
sition is completed. 



Whether the above phenomena about AD transition can be 
found for the general dynamics, say, for example, chaotic 
oscillators? To answer this, we have checked the transi- 
tion of the same network but for chaotic Rossler oscillators. 
The dynamics of a single oscillator is described by Fi (x) = 
[-Wiyi - Zi, uj 1 x 1 + 0.165?/;, Zi(Xi - 10) + 0.2], where to l is 
the natural frequency of the zth oscillator. In simulation, u>i 
is chosen randomly from the range [1,3]. The coupling func- 
tion is H(x) — x. Fig. 4 is a collection of the amplitude 
distributions observed during the transition, which repeat the 
phenomena of Landau-Stuart oscillators. Besides the specific 
network in Fig. 1, we have also tested the transition to CAD in 
other networks, including changing the network size and the 
degree distribution, adopting the clustered networks, and con- 
sidering the degree assortativity. The general finding is that, 
given the degree distribution is heterogeneous, the phenomena 
of staged transition and stair structures will be presented. 

To summarize, we have studied the transition to CAD in 
scale-free networks of nonlinear oscillators, and found the im- 
portant role of node degree played in AD generation. Since 
many practical systems where AD is importantly concerned 
possess heterogeneous degree distribution, our findings might 
give some new thoughts to the relevant studies, say, for exam- 
ple, the stability and evolution of ecological networks lfT8i[l9ll . 
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